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ABSTRACT 

We calculate all point symmetries of the Fokker - Planck equation in one-dimen- 
sional Euclidean space. General expression of symmetry group action on arbitrary solu- 
tion of Fokker - Planck equation is presented. We propose new notation for the group- 
theoretic analysis of PDE. The Lie prolongation formula is derived as an example of the 
new notation. 



1. The symmetries of the Fokker - Planck equation in one dimension 

The object of our considerations is a special case of Fokker - Planck equation, which describes evolu- 
tion of ID continuum of non-interacting particles imbedded in a dense medium without outer forces. The 
interaction between particles and medium causes combined diffusion in physical space and velocities space. 
The only force, which acts on particles, is damping force proportional to velocity. 

The 3D variant of this equation was investigated in our work [1]. In this work fiindamental solution 
of 3D equation was obtained by means of Fourier transform. 

In present work we choose another object of research. Our goal is to find aU symmetries of the 
Fokker - Planck equation. For this purpose we use the classic techniques of group-theoretic analysis of 
PDE, developed in fundamental works [2-5]. 

The Fokker - Planck equation in one dimension is 

dn dn dn , d^n 

where 

n - n{t, X, v) - density; 
t - time variable; 
X - space coordinate; 
V - velocity; 

a - coefficient of damping; 

k - coefficient of diffusion. 

The list of symmetries of the Fokker - Planck equation in one dimension follows. The calculations of 
symmetries are rather awkward. They are carried out to APPENDIX 2. 

Instead of classic (Z*" notation we use another ("S") notation. This notation is presented in AP- 
PENDIX 1. The Lie prolongation formula is derived as an example of the new notation. 
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The symmetries of the Fokker - Planck equation in one dimension are: 

. 3 



dn 



(2) 



where A is arbitrary solution of the (1) equation. 
Scaling of density 



V2- n 



dn 



The reason of symmetries (2-3) is Unearity of PDE (1). 
Time shift 



V3 



Space translation 



(3) 



(4) 



V4 = 



dx 



Extended Galilean transformation, which besides time and space coordinates affects the density 

dv dx 2k ^ dn 



(5) 



(6) 



Negative exponent transformation - it affects time and space coordinates, contains time-dependent common 
multiplier (negative exponent). It does not affect density. This transformation is evident from [1]. 



V6-e-|-a — + — 
dv ox 



(7) 



Positive exponent transformation - it affects time, space and density, contains time-dependent common 
multiplier (positive exponent). 



^ dv dx k dn 



(8) 



One-parameter groups, generated by vector fields Vi - V7, are enumerated in the following list. The 
list contains images of the point («, t, x, v) by transformation exp{£\i) 



Gl: 

G2: 
G3: 
G4: 



(n + eA, t, X, v); 
{e^n,t,x,v); 
in,t + £, X, v); 
(n, t,x + e, v); 



(9) 



(exp 



^ l^eiax + v) + - e\at + 1) 



n,t,x + et, V + e); 



(n, t,x + ee , v - ea e "'); 

1 



{n exp 



a 
~k 



£v+ - £^a e"' 
2 



,t,x + ee"', v + ea e"'). 
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Relatively nontrivial integration for G5 and G7 is carried out to APPENDIX 3 and APPENDIX 4. 

The fact, that G, are symmetries of PDE (1) means, that if /(/, x, v) is arbitrary solution of (1), the 
functions 



,,(1). 



,,(2). 



.,(3). 



,(4). 



f(t, X, v) + eA(t, X, v); 
e'f{t,x, v); 
/(/ - £, X, v); 
fit,x-£, v); 



(10) 



.,(5). 



exp 



eiax + v)-^e'^iat+l) 



f(t,x — £t,v — ey. 



.,(6). 



.,(7). 



f(t,x-ee-'",v + ea e'"'); 

1 



exp 



■e""| £v — e^a e"' 
k [ 2 



f{t,x-ee'",v-eae'"). 



where e - arbitrary real number , also are solutions of (1). Here A is another arbitrary solution of (1). 

We systematically replaced "old coordinates" by their expressions through "new coordinates". Note, 
that due to these replacements terms with in m^^' and m^^^ change their signs. 



We have trivial solution n - e°* at our disposal. If we act on this solution by transformations (10), we 
obtain 2 new solutions: 



n = exp 



a ( 1 , 

a/ - — £{ax + v) - - e {at + 1) 



(11) 



n — exp 



at e"'! £v e^a e" 

k [ 2 



General expression is 



u — e^^ exp 



a ( 1 2 

— f 5(ax + v) - - £5 (a? + 1) 



exp 



1 



2^ „at 



e"'\£yiv-e5)--£^a e 



(12) 



xf{t - £3, X- £4- £^t - £^e - £-je , V - £5 + £^a e - £-ja e ) + £iA{t, x, v). 
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DISCUSSION 

Looking at the list of all point symmetries of the Fokker - Planck equation in one-dimensional Euclidean 
space, we see, that there is no simple way to get, for example, fundamental solution of PDE, using these 
symmetries. We have not at our disposal such an instrument, as scaling of independent variables x, u, v. 

Indirect way of use of Galilean transformation (6) was demonstrated in [1]. The transformation was 
used for generalization of solution, which was obtained in the form of exponent of quadratic form of space 
coordinates and velocities with time dependent coefficients. 

There is need of further investigations of Fokker - Planck equation and its set of symmetries, which 
may lead to another physically interesting results. 
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APPENDIX 1 

We propose following notation. 

u" - the set of dependent variables; 

x' - the set of independent variables; 

9m" a . . . . 
^r-^ - u ,1 - first derivatives; 
ax' 

^ t . = u" ,ii - second derivatives. 
dx'dxi ■' 

The space of dependent and independent variables undergo small point transformation. The small 
displacements are: 

Su", Sx' - variations of independent and dependent variables. We call them "the components of infinitesi- 
mal transformation vector field" also. 

Variations Su", Sx^ are functions of independent and dependent variables. Therefore derivatives of 
variations on dependent and independent variables are: 

= iSu"),i , ^.f " = iSu"),H - derivatives of variations of dependent variables; 

dx' duP ^ 

dSx' ,. dSx' 

(Sx ),: , = (Sx ),n - derivatives of variations of independent variables. 



dxJ ' ■' ' duP 

In the same manner we write expressions for derivatives of the second and higher orders. 

Let us derive the Lie prolongation formula for variations of the first order derivatives. 

Due to point transformation of space of independent and dependent variables first derivatives undergo 
some transformation. To find this transformation, we write relation between differentials of independent 
and dependent variables 

du" = ^ dx\ (Al-1) 

dx' 

The same relation take place for variables after the transformation. Therefore we can take variation of 
(Al-1) and get 

3m" . 3m" 

S{du") = S(^)dx' + S(dx'). (Al-2) 
dx' dx' 

or 

3m" du" 
S(—)dx' S{du") - — Sidx''). (Al-3) 
dx' dx'' 

Variations of differentials are equal to 

Sidu") = d{Su") = Di{Su")dx'; (Al-4) 

S{dx'') = d{Sx'') = Di{Sx'')dx'; (Al-5) 

where symbol D, denotes the full partial derivative. 

3 3m" 3 

A = ^ + ^^ + ---. (Al-6) 



dx' dx' 3m" 
After substitution of (Al-4) and (Al-5) to (Al-3) we get 



3m" „ 3m" , 

S{^)^Di{Su'')-^Di{Sx''). (Al-7) 
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This is the famous Lie formula for the first derivatives (see ref. [1]). It is named "prolongation for- 
mula", because it prolongs the action of infinitesimal transformation to first derivatives. 

We shall use this formula in expanded form 

In the same way we get Lie prolongation formula for variations of the second order derivatives. 
We begin with 

^(^^) = D„iSiu"„ )) - D„iSxP). (Al-9) 

which is (Al-7) equation for u",^. 
or 

Siu",,^ ) - D„{S{u"„ )) - u"„p D„{5xP), (Al-10) 

where 5{u",k ) we take from (7). 

We expand expression for D„{S(u",j^ )) and obtain 

Siu ,1^ ) = uP ^ + D^({5x ),j ) + D^aSx ),p ) + (Al-11) 

From the (Al-7) we get 

a(<5(M",^)) dSu" dSx" dSx" ^ „ dSx"" 



and it follows 

9(^(m",^ )) ^ _ 9<^m" ^ _ a '^SxP dSxP y _ p a ^^X" 



In the similar way we have 

a,,^ , = Sjk, -> .. . D„iiSx'),j ) - D„ii5x'),k ). (Al-14) 

d(S(u",k)) d(S(u"',k)) « 

. ' , ' = => . ' . ) = ). (Al-16) 

d{{6uP),j) d{{5uP),j) 

„ = => ' ' . Z);„((^M^),r ) = D^iiSu"),y ). (Al-17) 

d{{SuP),y) d{{SuP),y) 

We collect aU these results and obtain Lie prolongation formula for variations of the second order de- 
rivatives: 
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5{li\^ ) = — U^kn. -U\pn, + ^ ) " (Al-18) 



dSx" 

) + i>„((^M'^),y ) - u"„p D„{sxn. 

We immediately use (Al-8), (Al-18) for our task. In our case the only dependent variable is u, inde- 
pendent variables are f, x, v. Expanding (Al-8), we have 

dn d d dn dn d d dn 

^ ^ = i^n) + ^ (Sn) - 5- (s- i^x) + ^ (Sx) ^) - (Al-19) 
ax ax an ax ax ax an ax 

dn d d dn dn d d dn 

- ^ (^ (^v) + ^ (^v) ^) - ^ (^ (St) + ^ (St) ^); 
dv ax an ax at ax an ax 

dn d d dn dn d d dn 

5^^^ {5n) + ^ (Sn) ^ - ^ (^ (Sx) + ^ (Sx) ^) - (Al-20) 
av av an av ax av an av 

dn d d dn dn d d dn 

- ^ (^ (Sv) + ^ (Sv) ^) - ^ (^ (St) + ^ (St) ^); 
av av an av at av an av 

dn d d dn dn d d dn 

^ 3- = ^ (Sn) + ^ (Sn) 37 - ^ (^ (Sx) + ^ (Sx) 3-) - (Al-21) 
at dt dn dt dx dt dn dt 

dn d d dn dn d d dn 

dv dt dn dt dt dt dn dt 

Expanding (Al-18), we have 

d^n d d^n d^n d d dn d^n d d dn 

^ ^TT = ^ (^«) - VlT ^JT ^^^^ + ^ ^^^^ JT^-JT^^JT ^^"^ + ^ ^"^^^ " (^1-22) 
dv^ an dv^ avax av an av av^ av an av 

d^n ^ d ^^^^ d ^^^^ dn ^ d^n dn ^ ^ ^ (S ) ^ {St + 

dtdv dv dn dv dv^ dx dn dv^ dv dn dv^ dt dn 

d^ d^ dn dn d^ 32 dn dn d^ 

+ ^ (Sn) + ^-3- (Sn) 3- + ^ (Sn) + ^ (Sn) ^)-^(^ (Sx) + 

av^ anav dv dv anav an^ dv ax dv^ 

d^ dn dn d^ 32 dn dn d^ d^ dn 

anav av av dv^ anav av at dv^ anav av 

dn dn ^ d^ (^x)-l- (Sx)^") ^" ^" ( iSv)+ (Sv)^") ^" ^" ( iSt) + 
dx dv dndv dri^ dv dv dv dndv dn^ dv dt dv dndv 

d^ dn d^n d d dn d^n d 

+ ^ (^0 ^)-^{^ iSx) + ^ {Sx) ^)-^{^ {Sv) + 
dn-^ dv dvdx dv dn dv dv^ dv 

d dn d^n d d dn 

+ ^ (^v) ^) - ^ (^ (^0 + Y- (^*) 3-)- 
dn dv dtdv dv dn dv 

These expressions are typical. It is obvious, that to manipulate such expressions by bare hands is very 
hard task. We use simple JavaScript program to generate these expressions. 
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APPENDIX2 



The infiinitesimal invariance criteria for equation (1) is 

dn dn dn dn dn d^n 

S{^) + ^ + v^(^) -aSv^- avS(^) - aSn - kd{^) = 0; (A2-1) 
at ax ax av av av^ 

dn 

We eliminate ^ from (Al) usmg original equation (1): 
at 

dn dn dn d^n 

^--v^ + av^ + an + k ^r— . (A2-2) 
dt dx dv dv^ 

Substitute (Al-19 - Al-22) to (A2-1), collect terms by different combinations of n derivatives and 
equate these terms to zero. We obtain following equations: 

dn dn 
dx 

-kv^St + k^Sx^Q (A2-3) 
dn^ dn^ 

dn dn 
dx dv 

d^ 3^ 3^ 

-kv .r-^ St - kv St + Ik Sx = Q (A2-4) 
dndv dndv dndv 

dn 
dx 

2 a d d d d^ 

—av — St + avn — St + av — 5x — an — Sx — kv ^-^ St + (A2-5) 
dv dn dv dn dv^ 

^ . d d d d 

+k^5x + v^ ^5t-v^Sx + v^St + Sv-^Sx = 
dv^ dx dx dt dt 



dn 
3v3 



dn d^n 
3v2 dv^ 



dn 

3^ 



dn d^n 
dv 



dn d^n 
dv dvdx 



d^ 32 

akv — St + k — Sv^O (A2-6) 
dn'^ dn'^ 



g2 

k'^^St = Q (A2-7) 
dn^ 



3^ g2 g2 g2 

lakv St + akn^St-k^Sn + 2k (A2-8) 

dndv dn^ dn-^ dndv 



2k^ St + 2k^Sv = (A2-9) 
dndv dn 
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an 



(A2-10) 



dn d^n 



dv dtdv 



dn 

a7 



8^ 



d^n 
dvdx 



^n 
dtdv 



2k — St = (A2-11) 
an 



3 3 3^ 3^ 3 3 

a^v^ ^St- a^vn — St + akv ^St + 2akn St - av^ — St + av — Sv- (A2-12) 
3v an 3v^ anav ax av 

3 3 32 32 3 3 

—av — St — an — Sv — aSv — 2k Sn + k —— Sv - v — Sv - — Sv-0 
at an anav av^ ax at 



3 3 , 32 3 3 3 

akv ^St-akn^5t + k^ ^St-kv^St + 2k^Sv-k^5t = (A2-13) 
av an av^ ax av at 



2k^Sx^0 (A2-14) 
av 



2k — St^0 
av 



(A2-15) 



a^vn ^ St- a^n^ ^ St + akn ^— r St - avn — St - av — Sn + 



dv 



dn 



3v2 



3x 



3v 



(A2-16) 



3 3 32 3 3 

+an ^ Sn — an TT- St — aSn — k zr-^ Sn + v zr- Sn+ zr- Sn — 0. 
dn dt av2 ax dt 



From (A2-10), (A2-11), (A2-14), (A2-15) we see, that Sx ^ Sxix,ty,St ^ StixJ). Using this, we 
simpUfy the rest of equations (A2-3 - A2-16). 

3 „ 3 _ 3 „ _ 3 
dx 



dx 



dt 



dt 
32 

A: ^ - 0; 
dn^ 



g2 g2 

-k^Sn + 2k .r-^ Sv = 0; 
dn^ dndv 



2k^Sv = 0; 
dn 



av ^ St + av ^ Sv — av ^ St — an ^ Sv 



dx 



dv 



dt 



dn 



(A2-17) 
(A2-18) 
(A2-19) 
(A2-20) 
(A2-21) 
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3^ d d 

—adv — 2k zr-zr- dn + k r 5v - V zr- 6v - zr- Sv = 0; 
onav dv^ dx dt 

d 3 3 

-kv^St + 2k^Sv-k^St^0; (A2-22) 
dx dv dt 

3 3 3 3 32 3 3 

—avn ^ St — av ^ Sn + an ^ Sn — an ^ St — aSn — k ^—r Sn + v ^ Sn + ^ Sn = 0. (A2-23) 
dx dv dn dt dv^ dx dt 

From (A2-20) we see, that Sv = Sv{x, v, t). Using this, we simplify the rest of equations (A2-17 - 
A2-23). 

2 3 3 3 3 

^St-v^Sx + v^St + Sv-^Sx = Q; (A2-24) 
dx dx dt dt 

32 

-k^Sn = 0; (A2-25) 
dn^ 

2 3 3 3 32 

—av ^ St + av zr- Sv - av zr- St - aSv - k ^ ^ Sn— (A2-26) 
dx dv dt dndv 

32 32 3 3 

—k ^ ^ Sn + k zr-^ Sv - V ^ Sv - zr- Sv = 0; 
dndv dv^ dx dt 

3 3 3 

-kv ^ St + 2k ^ Sv - k ^ St = 0; (A2-27) 
dx dv dt 

3 3 3 3 32 3 3 

-avn ^ St - av T— Sn + an z— Sn - an St - aSn — k ^r-^ Sn + v z^- Sn + ^ Sn = 0. (A2-28) 
dx dv dn dt dv^ dx dt 

From (A25) we conclude, that 

Sn:^A + nB; (A2-29) 

where A = A{x, v, t), B - B{x, v, t). Using this, we simplify the rest of equations (A2-24 - A2-28). 

, 3 3 3 3 

v^^St~v^Sx + v^St + Sv-^Sx = Q; (A2-30) 
dx dx dt dt 

, 3 3 3 32 3B 3 3 

-av^ ^ St + av ^ Sv - av ^ St - aSv + k Sv -2k ^ -v Sv - — Sv = 0; (A2-31) 
dx dv dt dv^ dv dx dt 

3 3 3 

-kv — St + 2k — Sv-k — St = 0; (A2-32) 
dx dv dt 

3 3fi 3 , d^B dB dB ^ 
-av^St-av^-a^St-k -I- v — + — = 0; (A2-33) 

dx dv dt dv^ dx dt 

dA d^A 3A 3A 

-av — -aA-k -— + v — + — = 0. (A2-34) 

dv dv^ dx dt 

Equation (A2-34) is simply Fokker - Planck equation for A. 
We solve (A2-30) and find Sv 

, 3 3 3 3 

Sv=-v^ — St + v — Sx-v — St+ — Sx. (A2-35) 
dx dx dt dt 

This gives for (A2-31 - A2-33) 
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dx 



St — av ^ 5t — a ^ dx — 2k 



dt 



dt 



dx 



St 



(A2-36) 



-lk — + v^—St + 2v^.^ 
dv ax^ atax 



9x2 



Sx + v z-^ St — 2v 



Sx- ^ Sx = '^\ 
dtdx dt^ 



3 3 3 

-5kv — St + 2k — Sx-3k — St = 0; 
ax dx at 



3 3B 3 d^B dB dB 

-av ^ St- av a ^St — k ^r-— -i- v — i- — = 0. 

dx dv dt dv^ dx dt 



(A2-37) 



(A2-38) 



Now we can collect similar terms by different degrees of v in (A2-37) and so split (A2-37) into two 
equations: 

3 



dx 



St = 0; 



3 3 
2k^Sx-2>k^5t = 0. 
dx dt 

From (A2-37a) we see, that St - St{t), which results in further simplifications 



3 . ^ . 
-av ^St-a^Sx-2k 
dt dt dv 



32 32 32 

+v — St - 2v z—- Sx- — Sx-0; 
3/2 3/3x 3?2 

3 3 
2kz^Sx-3kz^St^0; 
dx dt 



dB 3 d^B dB dB 

-av a — St — k + v ^ — I- ^ = 0; 



3v 3/ 



3v2 



3/ 



We integrate (A2-40) and find 



„ 3x 3 „ 
Sx^C + — z^St; 
2 dt 



where C - C{t). Further, we substitute this expression in (A2-39), (A2-41) and obtain 

3ax 32 3 3C , 35 3x d^ 3^ d^C 

- — - St - av — St - a — - 2k — — St-2v — St - — r- = 0; 

2 3/2 3/ 3/ 3v 2 3/3 3/2 3/2 

dB 3 d^B dB dB 

—av — a — St - k r- + v — I- rr— = 0. 

dv dt dv2 dx dt 

We integrate (A2-43) and obtain following expression (D = D(x, /)): 



B = D- 



1 f 3avx 3 



2k 



dC 3vx 3^ 



d^C 



^ ^^■St+ — —St + va — + — — St + v^z-rSt + v^. 
2 3f2 2 3f 3/2 3/3 3/2 3/2 



The substitution of (A2-45) to (A2-44), collecting and equating to zero terms by v gives 

l/2a^k- 



3 33 
3/ 3/3 



, 1 3^ 9 , 3C 1 3^ , 3^C dD 

3/4a^k-^x .^St+ l/2a^r> ^ - 3/4r^x ^ ^/ - l/2r' ^ + ^ = 0; 
3/2 3/ df^ 3/3 dx 



(A2-37a) 
(A2-37b) 

(A2-39) 



(A2-40) 
(A2-41) 

(A2-42) 

(A2-43) 
(A2-44) 

(A2-45) 

(A2-46) 
(A2-47) 
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3 3^ dD 
-l/2a ^St + ^St + ^=0. (A2-48) 



We integrate (A2-47) and obtain following expression {E - E{t)): 

3a^x^ 3^ xa^ 3C x d^C 



-Ilia — dt + ^dt + — =Q; (A2-52) 



^ = ^i^3^^^" 2^37-^3^^^-2137^] 
The substitution of (A2-49) to (A2-48), collecting and equating to zero terms by x gives 

-3/8a^r^ ^6t + 3/8r ' ^ - 0; (A2-50) 
3?3 3/5 

3^C 3'^C 

-l/2a2ri — + l/2fe-' ^ = 0; (A2-51) 

3 32 3£ 

— St+z-:dt+ — 

dt 3/2 3/ 
From (A2-46) and (A2-50) we conclude, that 

St ^ const ^Ci. (A2-53) 

From (A2-51) we conclude, that 

C ^C2 + C3t + C4e'" + €56-"'. (A2-54) 
From (A2-52) and (A2-53) we see, that 

E = Cg. (A2-55) 

We substitute expressions (A2-53 - A2-55) and obtain final expressions for variations (see (A2-35), 
(A2-42)): 

a^xn „ a^vn ^ „, avn ^ ^ . 

- - — C3 - — C^e"' - ^ <^3 + «Q + A; (A2-56) 

Sx^t C3 + C2 + C4e'" + €56-"'; (A2-57) 

^v = a C4 e^'-a C5 g-^' + Ca; (A2-58) 

St^Cu (A2-59) 

and these are final expressions for variations. 
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APPENDIX3 



In this appendix we integrate equations for Galilean transfomiation. 
The vector field of Galilean transformation is 

dv dx 2k ^ ^ dn' 
Characteristic differential equations are 

2k dn dt dv dx 

de ■■ 



We get first integrals 



an ax + v 1 t ' 



t-tQ, V - S + Vq, X — £t + Xq. 



The only nontrivial equation is 



(ax + v) ds ■ 



2k dn 



a n 



{a{£t + Xq) + £ + Vq) de — — ■ 



2k dn 



a n 



We integrate (A3-5) and get 



1 2 2k n 
— £ {at + \) + £{axQ -I- Vq) = ln{ — ); 



Finally, we get following expression for n 

a (I 

n^noexp -- - 



a no 



£ {at+l) + £{axo + Vq) 



(A3-1) 

(A3-2) 

(A3-3) 

(A3-4) 
(A3-5) 

(A3-6) 
(A3-7) 
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APPENDIX4 



In this appendix we integrate equations for positive exponent transformation. 
The vector field of positive exponent transformation is 



^ ^ dv dx k dn f 



Characteristic differential equations are 



de — —e 



k dn dt _„, dv 
a^v n 



We get first integrals 
The only nontrivial equation is 



t — tQ-, v-Vo + eae"'; x-Xo + ee"'. 



-e'"v de ■■ 



k dn 



-e'"(vo + ea e"') de ^ 



k dn 



We integrate (A4-5) and get 



-e^'ievo+l-e^a e"') = - ln(-y, 
2 



k n 
no 



Finally, we get following expression for n 



n- no exp 



^ e-\£Vo+-e^a e"' 



(A4-1) 

(A4-2) 

(A4-3) 

(A4-4) 
(A4-5) 

(A4-6) 
(A4-7) 



